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Abstract 

We perform an investigation of the MSSM parameter space at the Fermi 
scale taking into account available accelerator, non-accelerator and cosmological 
{vq . constraints. Extra assumptions about upper bounds for some of the SUSY 

' particles are also imposed. We show that a non-observation of the SUSY dark 

matter candidates with a high-accuracy dark matter detector, such as 73 Ge, 
under above-mentioned assumptions can exclude large domains of the MSSM 
parameter space and, for instance, can make especially desirable collider search 
for light SUSY charged Higgs boson. 

(N ; 

! As already well-known [jl], [3[], a direct dark matter search for neutralinos, light- 

est SUSY particles (LSP), is complimentary to high energy searches for SUSY with 
colliders. The direct dark matter search is able to give information which is not 
available from collider physics. 
Q\ . It has been found, that future-possible SUSY spectrum restrictions from colliders, 

especially such as upper mass bounds for some of the SUSY particles, can strongly 
enforce the importance of the direct dark matter detection experiments. The main 
point is that this hypothetical information leads to some important restrictions of the 
expected event rates for direct dark matter detection. In some cases the restriction 
is a clear lower limit for this rate. Another obvious example of such restriction is 
j> ! an upper bound for the LSP, which appears only due to the fact that all sfermions 

and gauginos are not lighter than the LSP. Thus, these high-energy bounds bring in 
restrictions on the parameter space which can result in a lower limit for the expected 
rate. Such effect is obtained mostly when the sfermion masses are restricted, but the 
effect appeared quite strongly especially in the Higgs sector of the MSSM ||. 

The question — to what extent the accelerator searches are able to affect future 
prospects for direct dark matter search provided the dark matter particle candidate is 
LSP — was investigated comprehensively in the minimal supergravity SUSY models 
@) ||) 0' H- I n this paper we explore the MSSM parameter space at the weak scale in 
the most phenomenological way, relaxing completely any constraints following from 
the unification assumption. On the other side we hold all restrictions from the age of 
the Universe, accelerator SUSY searches, rare FCNC b — > 57 decay, etc |7|, |8|, [£|. 

The MSSM parameter space in our approach is determined by entries of the mass 
matrices of neutralinos, charginos, Higgs bosons, sleptons and squarks. To specify 
the parameter space we give all relevant mass matrices below. The one-generation 

laboratory of Nuclear Problems, Joint Institute for Nuclear Research, Moscow region, 141980 
Dubna, Russia. E-mail: bedny@nusun.jinr.ru 



1 



squark and slepton mass matrices (using the notation of the third family) are given 
by I 
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where Syy = sin 2 W and tan/3 = (H®)/^®}. In the basis, the chargino 

mass matrix is 

_ / M 2 V2«V sin ^ ' 

Two unitary 2x2 matrices U and V are required to diagonalize the chargino mass- 
squared matrix = VX^XV^ 1 = U*XX^(U*)" 1 . The two mass eigenstates are 

denoted by xf an d xt- m the B-W^-H^-H® basis, the neutralino Majorana mass 
matrix is 
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where s@ = sin /3, cp = cos /3, etc. A 4 x 4 unitary matrix Z is required to diagonalize 
the neutralino mass matrix Ai~o = Z*YZ~ 1 where the diagonal elements of M.~< 

X. X. 

be either positive or negative. The CP-even Higgs mass matrix has the form [y 
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Here = cos 2 /3, = sin 2 /3 and </ s /(m 2 ,m|) = 2— m ^ m | ln^. Neutral CP-even 
Higgs eigenvalues are 
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i{# u + # 22 ± ^/(tf n + # 22 ) 2 - 4(H U H 22 - H^j}, 



The mixing angle a H is obtained from 
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The mass of the charged Higgs boson is given by m^ H = m 2 ^ + M\ + cuA c h. 

Therefore free parameters are: tan j3 is the ratio of neutral Higgs boson vacuum 
expectation values, \i is the bilinear Higgs parameter of the superpotential, Mi, M 2 
are soft gaugino masses, Ma is the CP-odd Higgs mass, m~, m~, m~ are squared 

squark mass parameters for the 1st and 2nd generation, m~, m~ are squared slepton 

L E 

mass parameters for the 1st and 2nd generation, m~ , m~, m- are squared 3rd 

Q3 ± B 

generation squark mass parameters, m~ , m~ are squared 3rd generation slepton 
mass parameters and A t , A , A T are soft trilinear couplings for the 3rd generation. 
With these parameters one completely determines the SUSY spectrum and the MSSM 
coupling constants at the Fermi scale. 

A dark matter (DM) event is elastic scattering of a relic DM neutralino from a 
target nucleus producing a nuclear recoil which can be detected by a suitable detector 
||. The differential event rate in respect to the recoil energy is the subject of exper- 
imental measurements. The rate depends on the distribution of the DM neutralinos 
in the solar vicinity and the cross section of neutralino-nucleus elastic scattering. 
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In our analysis we use the so-called total event rate R which is integrated over 
recoil energies, and which is useful for searching for domains with extreme rates. 



We follow our paper |TT|, where on can find all relevant formulas and astrophysical 
parameters. To calculate the event rate we use for the relic neutralino mass density 
and for the escape neutralino velocity commonly accepted values 0.3 GeV/cm 3 and 
600 km/s, respectively. Their experimental variations can slightly change R but leave 
the dependence of R on the MSSM parameters (Fig. [I] and||) unaffected. To compare 
our results with sensitivities of different dark matter experiments we calculate also 
the total cross section for relic neutralino elastic scattering on the nucleon. 

The present lifetime of the Universe implies an upper limit on the expansion rate 



and correspondingly on the total relic abundance. One finds |T2[ that the contribution 
of each relic particle species x has to obey ^l x h^ < 1, where h is the Hubble constant 
and the relic density parameter Q x = p x /p c is the ratio of the relic neutralino mass 
density p x to the critical one p c = 1.88-10 _29 hQg-cm -3 . Assuming that the neutralinos 
form a dominant part of the dark matter in the Universe one obtains a lower limit on 
the neutralino relic density. We restrict our analysis with the cosmological constraint 
0.025 < fi x /ip < 1 10) 0; O- We calculate £l x h^ following the standard procedure on 



the basis of the approximate formula [13|, |14 |. We take into account all channels of 



the x ~ X annihilation. Since the neutralinos are mixtures of gauginos and higgsinos, 
the annihilation can occur both, via s-channel exchange of the Z° and Higgs bosons 
and t-channel exchange of a scalar particle, like a selectron [15|. This constrains the 
parameter space || [L3|, [1^ . 



Another stringent constraint is imposed by the branching ratio of b — > sj decay, 
measured by the CLEO collaboration to be 1.0 x 10 -4 < B(b — > 57) < 4.2 x 10~ 4 . 
In the MSSM this flavor changing neutral current process receives contributions from 
H —t, x —t and g-q loops in addition to the standard model W— t loop. These also 
strongly restrict the parameter space [17] . 



The masses of the supersymmetric particles are constrained by the results from the 
high energy colliders LEP at CERN and Tevatron at Fermilab. This imposes relevant 
constraints on the parameter space of the MSSM. We use the following experimental 



restrictions for the SUSY particle spectrum in the MSSM 181: M~+ > 65 GeV for 

A 2 



the light chargino, M-+ > 99 GeV for the heavy chargino, M^o^^ > 45, 76, 127 GeV 
for non-LSP neutralinos, respectively; M v > 43 GeV for sneutrinos, M^ R > 70 GeV 
for selectrons, Mq > 210 GeV for squarks, > 85 GeV for light top-squark, 
M H o > 79 GeV for neutral Higgs bosons, M C h > 70 GeV for charged Higgs boson. 

In our numerical analysis the parameters of the MSSM are randomly varied at the 
Fermi scale in the intervals given below 
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Figure 1: Total event rate in 73 Ge versus MSSM parameters M 2 , tan/3, /i, M A) m 2 Q l2 
(labeled as "sq2") and m\ (as "sl2"). 

For each parameter some number, between and 1, was defined by the random 
generator from the CERN Library. With this number and lower and upper bounds 
given above the random value of the MSSM parameter was calculated by means of 
linear function. We stopped running when the lower border of the rate R stopped 
move down and only density of points continued to increase. In our main scan we 
have tested about 10 8 models, and only about 10 5 models passed all constraints. 
For simplicity, for other sfermion mass parameters we used the relations m~ = 

m ~ = m ~ , m ~ = m 2 r , m~ = m ~ = m%„ , m ~ = mh o . The parameters A\, and 
A T are fixed to be zero. We consider the domain of the MSSM parameter space, in 
which we perform our scans, as quite spread and natural. Any extra expansion of it 
like, for example, using —10 TeV< M 2 < 10 TeV, etc, of course, can be possible, but 
should be considered as quite unnatural in the framework of the idea of SUSY. 

Some of results of the main (without extra constraints) scan are presented in 
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Light Higgs mass in GeV Charged Higgs mass in GeV 

Figure 2: Total event rate in 73 Ge versus masses of squarks, sleptons, light chargino 
(practically the same scatter plot one obtains for next-to-lsp neutralino), heavy 
chargino, light CP-even neutral Higgs boson and charged Higgs boson (the same 
is for heavy Higgs boson). 

Figs. [l]-f3]as scatter plots. The main feature we paying our attention to is the presence 
of a lower bound for the total event rate R. An absolute minimum value of about 10~ 6 
events/day/kg in a 73 Ge detector is obtained in the above-mentioned domain of the 
MSSM parameter space. There is a clear growth (up to one order of magnitude) of 
the lower bound only with tan/? (Fig. |IJ). In all other cases there is a decrease of the 
lower bound, the decrease is most sharp with \fi\, M A (about 5 orders of magnitude) 
and Mq i 2 (Fig. [1]) and with squark mass M~, heavy chargino mass and charged 
Higgs boson mass Men (Fig. 0). 

Figure |3] (upper panel) gives expectations for the total event rate R obtained by 
means of scanning the parameter space at the Fermi scale without any extra limita- 
tions on the SUSY particle spectrum. There is a lower bound, which decreases with 
mass of the LSP and reaches an absolute minimum of about 10~ 6 events/day/kg in a 
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Figure 3: Total event rate in 73 Ge in events/day/kg (upper panel) and cross section 
of scalar elastic scattering of LSP (WIMP) off neutron in GeV -2 (lower panel) versus 
mass of LSP (in GeV). 



73 Ge detector in the region of LSP masses 600-700 GeV. The phenomeno logically al- 
lowed masses of the LSP are spread from about 5 GeV till about 800 GeV. Considering 
both scatter plots in Fig. § one also can conclude that, due to the practical absence 
of a lower bound for the scalar cross section of WIMP-neutron interaction, the lower 
bound for the rate is mainly established by the spin- dependent interaction, which in 
contrast to the scalar interaction is associated with an about 3-order-of-magnitude 
larger lower bound for WIMP-nucleon cross section. 

The existence of the lower bound for the event rate itself and the variation of 
the bound with the MSSM parameters and masses of the SUSY particles allow us 
to consider prospects to search for dark matter under special assumptions about 
restricted values for the above-mentioned parameters and masses. To this end we 
have performed a number of extra scans taking into account extra limitations on 
single squark mass (M sq < 250, 230 GeV), light neutral CP-even Higgs boson mass 
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Figure 4: Different lower bounds for the total event rate in 73 Ge (events/day /kg) 
versus mass of the LSP (GeV). Here M CH H1 denote masses of the squark, the 
charged Higgs boson and the light neutral CP-even Higgs boson respectively. Heavy 
chargino mass is denoted as M ch _ os . "Full" corresponds to the lower bound obtained 
from main (unconstrained) scan, and " Light spectrum" denotes the lower bound for R, 
which is obtained with all sfermion masses lighter then about 300 GeV. The horizontal 
dotted line represents expected sensitivity for the direct dark matter detection with 
GENIUS. 



(M H1 < 80, 100, 120 GeV), charged Higgs boson mass (M CH < 150, 200 GeV) 
and heavy chargino mass (M ch _ os < 250 GeV). We also considered the situation 
where masses of all superpartners not exceed 300-400 GeV. All corresponding curves 
together with the absolute lower bound from the unconstrained scan are depicted in 
Fig. |. 

A restriction of the single (light) squark mass to be quite small (M sq < 230 GeV) 
as well as another assumption that all sferminos masses not exceed 300-400 GeV, 
put upper limits on the mass of the LSP and therefore do not permit R to drop very 
deeply with increasing LSP mass. Furthermore in both cases the lower bound for the 
rate is established for all allowed masses of the LSP at a level of 10~ 3 events/kg/day. 
This value for the event rate we consider as an optimistic sensitivity expectation for 



high-accuracy future detectors of dark matter, such as GENIUS [[19], g(J . Practically 
the same lower bound one obtains under the assumption that both charginos are quite 
light (M ch _ os < 250 GeV). In this case LSP masses do not exceed 250 GeV and lower 
rate bound is 2 5 x 10 -3 events/kg/ day. 

One can see that the mass of the light neutral CP-even Higgs boson M H1 , perhaps 
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Figure 5: WIMP-nucleon cross section limits in pb for scalar interactions as function 
of the WIMP mass in GeV. Filled circles present our calculations with light SUSY 
spectrum. Filled triangles give the cross section with assumption of M CH < 200 GeV. 

most easily measurable experimentally due to its smallness, has unfortunately only a 
very poor restrictive potential. Already for the mass value M H1 < 80 GeV, which is 
practically equal to the experimental border (M^ p = 79 GeV), the lower bound for 
the event rate in a detector of 73 Ge is far below 10~ 3 events/kg/day. The situation 
looks most promising when one limits the mass of the charged Higgs boson. 

From Fig. ^| one can conclude that in the charged Higgs boson low-mass region 
other masses of both CP-even and CP-odd Higgs bosons are also restricted from 
above. Therefore coupling constants of the scalar neutralino-quark interaction, which 
contain terms with ^ — factors, are not suppressed enough and the rate can not 

decrease significantly. The lower bound of the rate increases when the mass M CH 
decreases (Fig. g) and for M CH < 200, 150 GeV reaches values of ~ 10 -2 , ~ 10 _1 
events/kg/day, respectively practically for all allowed masses of the LSP. These val- 
ues can be reached not only with GENIUS (10 -2 events/kg/day), but also with some 
other near- future direct dark matter detectors The fact also is confirmed by 

Fig. [5], where cross section limits for the WIMP-nucleon scalar interactions for differ- 
ent experiments are presented as functions of the WIMP mass. Light-filled circles in 
Fig. |5] give the scalar cross sections, calculated under the assumption that the SUSY 
spectrum is quite light (see Fig[|). Filled triangles give the cross section, obtained 
with charged Higgs boson mass restriction M CH < 200 GeV. 

Therefore if it happened, for instance, that either the SUSY spectrum is indeed 
light, or the charged Higgs boson mass indeed does not exceed 200 GeV, in both cases 
at least the GENIUS experiment should detect a dark matter signal. If we consider a 
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Figure 6: WIMP-nucleon cross section limits in pb for scalar interactions as function 
of the WIMP mass in GeV. Filled circles present our calculations with light SUSY 
spectrum and tan/3 > 20. Filled triangles give the same as filled circles, but for 
tan ft > 40. 

more complicated condition, for example, assuming that the SUSY spectrum is quite 
light and simultaneously that tan/3 is quite large, then not only GENIUS, but also 
CDMS and HDMS ^TJ will possess very good prospects to detect a dark matter 
signal. This situation is illustrated in Fig. || where besides cross section limits for the 
WIMP-nucleon scalar interactions for different experiments are given calculations for 
the case of a light SUSY spectrum with extra assumptions tan ft > 20 (filled circles) 
and tan ft > 40 (filled triangles). 

Therefore the obtained correlations between lower limit for the event rate R and 
some masses of SUSY particles give good prospects for direct dark matter detection 
with next-generation detectors. The prospects could be very promising if from collider 
searches one would be able to restrict the mass of the charged Higgs boson at a level of 
about 200 GeV (light Higgs sector). The observation, due to its importance for dark 
matter detection, could serve as a source for extra efforts in searching for charged 
Higgs boson with colliders. Considered together, both these experiments, collider 
search for charged Higgs and dark matter search for SUSY LSP, become very decisive 
for a verification of SUSY models. 

Otherwise complete non-observation of any dark matter signal with very sensitive 
dark matter detectors in accordance with Fig. would exclude, for example, a 
SUSY spectrum with masses lighter then 300-400 GeV as well as light SUSY spectrum 
with large tan ft (Fig. |6]), charginos with masses smaller then 250 GeV (Fig. |j), charged 
Higgs boson with M CH < 200 GeV, and therefore any possibility for the entire light 
Higgs sector in the MSSM (Fig. g). 

The last case is in particular interesting, because if the light charged Higgs boson 
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is excluded by GENIUS, then either it will be rather unpromising to search it for 
with colliders, or any positive result of a collider search brings strong contradictions 
in the MSSM approach to dark matter detections and/or collider SUSY searches. 

Finally we would like to comment the following. Unfortunately the MSSM param- 
eter space is huge and to obtain some reliable feeling, concerning, for example, the 
expected rate of dark matter detection when all relevant experimental and cosmo log- 
ical constraints are taken into account, one has nothing but this statistical numerical 
method 0, g, §, 0, |I6|. 



This method allows lower and upper bounds for any observable to be estimated, 
and to make conclusions about the prospects for dark matter detection with modern 
or near-future high-accuracy dark matter detectors. The larger the amount of points 
which confirms such a conclusion the better. The conclusions we made here are based 
on hundreds of thousand of points which passed all constraints. Of course, we have 
no proved protection against peculiar choices of parameters which could lead to some 
cancellation and to small cross sections even if Higgs masses are small. Nevertheless, 
the probability of these choices is very small (about 1/100000), otherwise we should 
already meet them with our random scanning. On the other side, if these peculiar 
choices exist and one-day would manifest themselves, this would be a very interesting 
puzzle, because it would be some kind of fine tuning of parameters, which requires 
strong further development of our understanding of the theory. 

The results of this paper may be considered as a good example of the complemen- 
tarity of modern accelerator and non-accelerator experiments looking for new physical 
phenomena. 

We thank S.G. Kovalenko and L. Baudis for helpful discussions. The investigation 
was supported in part (V.A.B.) by Grant GNTP 215 from Russian Ministry of Science. 
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